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Abstract. Let d > 3. We consider the global Cauchy problem 
for the generalised Navier-Stokes system 

dtu + {u ■ V)w = -D'^u - Vp 

V • u = 

m(0, x) = Uq{x) 

for u : R+ X R"^ ^ R"* and p : R+ X R"* ^ R, where uq-.R'^ ^ R'^ 
is smooth and divergence free, and D is a Fourier multiplier whose 
symbol m : R'^ R+ is non-negative; the case m(^) — j^j is 
essentially Navier-Stokes. It is folklore (see e.g. [S]) that one has 
global regularity in the critical and subcritical hyperdissipation 
regimes to(C) = |^|" for a > We improve this shghtly by 

establishing global regularity under the slightly weaker condition 
that to(^) > for all sufficiently large ^ and some 

non-decreasing function g : R+ — > R+ such that sgfly ~ 
In particular, the results apply for the logarithmically supercritical 
dissipation m(0 := ICl"^/ log(2 + |e|2)i/4. 



1. Introduction 

Let d > 3. This note is concerned with solutions to the generahsed 
Navier-Stokes system 

dtU + [u ■ V)u = —D^u — Vp 

V ■ M = (1) 

u(0, x) = Uo{x) 

where u : R"*" x R'^ —>■ TH^, p : R"*" x R'^ R are smooth, and Uq : 
R'^ — > R"' is smooth, compactly supported, and divergence-free, and D 
is a Fourier multiplieiEl whose symbol m : TV^ —>■ R"*" is non-negative; 
the case m(^) = |^| is essentially the Navier-Stokes system, while the 
case m = is the Euler system. 
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"'^Thc exact definition of the Fourier transform is not relevant for this discussion, 
but for sake of concreteness one can take /(^) := /(a^)e~"^ dx. 
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For d > 3, the global regularity of the Navier- Stokes system is of course 
a notoriously difficult unsolved problem, due in large part to the super- 
critical nature of the equation with respect to the energy E{u{t)) := 
/ad dx. This supercriticality can be avoided by strengthening 

the dissipative symbol ?7t,(^), for instance setting m(^) := |^|" for some 
a > 1. This hyper- dissipative variant of the Navier-Stokes equation 
becomes subcritical for a > (and critical for a = ^j^), and it is 
known that global regularity can be recovered in these cases; see [5] for 
further discussion. For 1 < a < only partial regularity results are 
known; see ^ for the a = 1 case and [5J for the a > 1 case. 

The purpose of this note is to extend the global regularity result very 
slightly into the supercritical regime: 

Theorem 1.1. Suppose that m obeys the lower bound 

m(0 > l^l^'^'^^Vgm (2) 

for all sufficiently large |^|, where g : is a non- decreasing 

function such that 

Then for any smooth, compactly supported initial data Uq, one has a 
global smooth solution to ([T]). 



Note that the hypotheses are for instance satisfied when 

^(^):=|^|(^+2)/4/i^gi/4^2 + |er), (4) 

thus 

W/2(2-A)' 

Analogous "barely supercritical" global regularity results were estab- 
lished for the nonlinear wave equation recently in [S], [H], 

The argument is in fact quite simple, being based on the classical energy 
method and Sobolev embedding. The basic point is that whereas in 
the critical and subcritical cases one can get an energy inequality of 
the form 

^tlkWII^HR'*) ^ C'a(t)||«(t)||^,(j^d) 
for some locally integrable function a{t) of time, a constant C, and 
some large k, which by Gronwall's inequality is sufficient to establish 
a suitable a priori bound, in the logarithmically supercritical case (jlj) 
one instead obtains the slightly weaker inequality 

< Ca{t)\\u{t)\\]jk^^d) log(2 + \\u{t)\\jjk^-^d)) 

(thanks to an endpoint version of Sobolev embedding, closely related 
to an inequality of Brezis and Wainger [2]), which is still sufficient to 
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obtain an a priori bound (though one which is now double-exponential 
rather than single-exponential, cf. [IJ). 

Remark 1.2. It may well be that the condition can be relaxed fur- 
ther by a more sophisticated argument. Indeed, the following heuristic 
suggests that one should be able to weaken ([3]) to 7^)^ ~ ^^^^ 
allowing one to increase the 1/4 exponent in (jlj) to 1/2. Consider a 
blowup scenario in which the solution blows up at some finite time T^,, 
and is concentrated on a ball of radius 1/N{t) for times < t < T*, 
where N{t) oo as t ^ T^:. As the energy of the fluid must stay 
bounded, we obtain the heuristic bound u{t) = 0{N{tY/'^) for times 
< t < T*. In particular, we expect the fluid to propagate at speeds 
0{N{tY/'^), leading to the heuristic ODE |^ = 0{N{tY/^) for the 
radius 1/N{t) of the fluid. Solving this ODE, we are led to a heuristic 

2 

upper bound N{t) = 0{{T^ — 1)'^+^) on the blowup rate. On the other 
hand, from the energy inequality 

2 I f \Du{t,x)\'^ dxdt < [ \uo{x)\'^ dx 

Jo Jr'' iR<* 

one is led to the heuristic bound 



dt < oo. 



]V(t)(rf+2)/2^(]V(t))2 

2 

This is incompatible with the upper bound N{t) = 0{{T^ — t)^) if 
Ii* sg{s)^ ~ Unfortunately the author was not able to make this 
argument precise, as there appear to be multiple and inequivalent ways 
to rigorously define an analogue of the "frequency scale" N(t), and all 
attempts of the author to equate different versions of these analogues 
lost one or more powers of g{s). 

To go beyond the barrier 7^7)2 ~ ^ (with the aim of getting closer 

to the Navier-Stokes regime, in which g(s) = s^^^ in three dimensions), 
the above heuristic analysis suggests that one would need to force the 
energy to not concentrate into small balls, but instead to exhibit tur- 
bulent behaviour. 

1.3. Acknowledgements. The author is supported by NSF Research 
Award DMS-0649473, the NSF Waterman award and a grant from the 
MacArthur Foundation. 



2. Proof of theorem 



We now prove the theorem. Let A; be a large integer (e.g. k := lOOd 
will suffice). 
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Standard energy method arguments (see e.g. [1]) show that if the 
initial data is smooth and compactly supported, then either a smooth 
H°° solution exists for all time, or there exists a smooth solution up 
to some blowup time < T^, < oo, and ||'u(t)||/^fc(pid) oo as t — > T*. 
Thus, to establish global regularity, it suffices to prove an a priori bound 
of the form 

||M(^)||Hfc{R'*) < C{k,d, \\uo\\H''iRd),T,g) 

for all < t < T < cx) and all smooth H°° solutions u : [0, T] x R'^ ^ 
R*^ to ([I]), where C{k,d, ||mo||h'={R'*); S') is a quantity depending only 
on k, d, ||iio||//fc(R,d), T, and g. 

We now fix uq, u, T, and let C denote any constant depending on k, d, 
||Mo||/i-fc(R,d), T, and g (whose value can vary from line to line). Mul- 
tiplying the Navier-Stokes equation by u and integrating by parts, we 
obtain the well-known energy identity 

dt / \u{t,x)\^ dx = -2a{t) 

where 

a{t) := \\Du\\l2(^^,^ (5) 

(note that the pressure term Vp disappears thanks to the divergence 
free condition); integrating this in time, we obtain the energy dissipa- 
tion bound 

[ a{t) dt < C. (6) 
Jo 



Now, we consider the higher energy 

Ek{t):=y^! \VH{t,x)\' dx. (7) 

Differentiating ([7]) in time and integrating by parts, we obtain 

dtEkit) = -2 V ||V^Z^n(t)||2 2(j^,)-2 V / V^u{t,xyV^{{u-W)u){t,x)dx- 

again, the pressure term disappears thanks to the divergence-free con- 
dition. For brevity we shall now drop explicit mention of the t and x 
variables. 

We apply the Leibniz rule to V^{{u ■ V)^). There is one term involving 
j + 1*'' derivatives of u, but the contribution of that term vanishes by 
integration by parts and the divergence free property. The remaining 
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terms give contributions of the form 

3=0 I<jl,j2<j-jl+j2=j + l 

where 0{'V^u'V^^u'V^^u) denotes some constant-coefficient trihnear com- 
bination of the components of V%, V^^u, and V^^u whose exphcit form 
is easily computed, but is not of importance to our argument. We can 
integrate by parts using D and and then use Cauchy-Schwarz to 
bound 



0{V^uV^'uV^^u) dx < ||(l+Z))V%|U2(Kd)||(l+D)-^(C(V^'iMV^'2M))||i2(Rd). 

By the arithmetic mean-geometric mean inequahty we then have 

f 0{V'uV''uV^'u) dx < c||(l+D)V%||i2(K.)+-||(l+/^)"'(0(V^'^nV^'^^i))|li2(R.) 

for any c > 0. Finally, from the triangle inequality, and the fact 
that D commutes with V-', we have 

11(1 + D)V%||i.(i,.) < Ci\\VWu\\l2^^,^ + E,). 

Putting this all together and choosing c small enough, we conclude that 

i<ii<i2<fc:ji+i2<fc+i 

(8) 

To estimate this expression, we introduce a parameter > 1 (depend- 
ing on t) to be optimised later, and divide {1 + D)~^ = {1 + D)~^P<n + 
(1 -f D)~^P^j^, where P<n and P^n are the Fourier projections to the 
regions : 1^1 < N} and : 1^1 > N}. 

We first deal with the low-frequency contribution to ([H])- From Plancherel's 
theorem and we see that 

\\{l+D)-'P^^{0{V^^uV^-u))h2(^^.)<Cg{N)\\{Vy^ 

where (V)"^"'^^-*/^ is the Fourier multiplier with symbol (|^)~(''+2)/4^ 
where (^) := (1 -|- Applying Sobolev embedding, we can bound 

the right-hand side by 

< (iV) 1 1 1 V^'^ n 1 1 V-^'^M 1 1 1 i4d/(3d+2) (Rd) . 

By Holder's inequality and the Gagliardo-Nirenberg inequality, we can 
bound this by 

< C(7(Ar)||VM|L4./,.+2)(R.)||V^'^+^'^-in|U2(K.) 

which by ([7]) is bounded by 
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Next, we partition 

II VM||^4d/(d+2)(p{_d) < \\'VP<NU\\l4d/(d+2)(^^d-) + II VP>ArM||^4d/(d+2)(p_d). 

From Sobolev embedding and Plancherel, ([2]) and we have 
||VP<^M|L4d/(d+2)(Kd) < C||(V)('^+')/^P<jvm|| 

<Cg{N){l + a{t)y/'. 
Meanwhile, from Sobolev embedding we have 

II VP>ArM||^4d/(d+2)(Rd) < jy-^fc 

(say) if k is large enough. Putting this all together, we see that the 
low-frequency contribution to ([H]) is 

< Cg{NfE,[g{N)\l + a{t)) + -^E,]. 

Next, we turn to the high-frequency contribution to From Plancherel, 
Holder's inequality, and ([7j) we have 

11(1 + D)-^P>jv(C(V^'^MV^'2M))|U2(Rd) < C(?(Ar)Ar-('^+2)/4|||yji^||y,.^|||^ 

while from Sobolev embedding and (^^) we see (for k large enough) that 

llV^'^nlUoc(Rd) < CE'J'. 

Thus the high-frequency contribution to ([8]) is 

< C(7(iV)2iV~(''+2)/2^2^ 

Putting this all together, we conclude that 

dtEk < Cg{NfEk[g{N)\l + a{t)) + ^E,]. 

We now optimise in N, setting := 1 + P^,, to obtain 

dtEk<Cg{l + E,)^Ek{l + a{t)). 

From ([6]), ([3]) and separation of variables we see that the ODE 

dtE = Cgil + EYE{l + a{t)) 

with initial data P(0) > does not blow up in time. Also, from ([7]) 
we have Pfc(O) < C. A standard ODE comparison (or continuity) 
argument then shows that Ek{t) < C{T) for all < t < T, and the 
claim follows. 

Remark 2.1. It should be clear to the experts that the domain R*^ here 
could be replaced by any other sufficiently smooth domain, e.g. the 
torus R°'/Z'^, using standard substitutes for the Littlewood-Paley type 
operators P<n, P>n (e.g. one could use spectral projections of the 
Laplacian). We omit the details. 
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